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ABSTRACT 

The discovery of Lya emission from galaxies at redshifts beyond z ~ 6.5 should not be naively interpreted 
as implying that the intergalactic medium (IGM) had been reionized at higher redshifts. We show that 
a cluster of faint undetected sources around each observed galaxy generates an HII region sufficiently 
large to allow transmission of the galaxy's Lya line prior to reionization. We also show that quasars may 
contribute a significant fraction of the ionizing photons to HII regions around galaxies with a velocity 
dispersion larger than ~ 100 km s~^. These contributing quasars are not usually seen due to the small 
fraction of time they spend in a luminous phase. 

Subject headings: cosmology: theory - galaxies: formation 
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1. INTRODUCTION 

The identification of a Gunn-Peterson trough in the 
spectra of the most distant quasars at redshifts of z '--^ 6.3- 
6.4 (Fan et al. 2004) hints that the reionization of cosmic 
hydrogen was completed only around that time, about a 
billion years after the big bang (White et al. 2003; Wyithe 
& Loeb 2004). On the other hand, the recent discovery of 
Lya emission from galaxies beyond a redshift of z ~ 6.5 
(Hu et al., 2002; Kodaira et al., 2003; Rhoads et al. 2004; 
Kurk et al. 2004; Stern et al. 2004) was naively interpreted 
as an indication that the surrounding intergalactic medium 
(IGM) was highly ionized, in apparent conflict with the 
quasar data. Although a galaxy is often surrounded by 
an HII region of ionized hydrogen (Haiman 2002; San- 
tos 2004), its Lya line could be heavily suppressed by the 
red damping wing of the resonant Lya absorption in the 
surrounding IGM (Loeb, Barkana & Hernquist 2004; Cen, 
Haiman & Messinger 2004; Ricotti et al. 2004). Cluster- 
ing of sources may result in HII regions that are larger 
than expected for a single galaxy, rendering the optical 
depth from the red damping wing of the IGM irrelevant 
(Gnedin & Prada 2004; Furlanetto, Hernquist & Zaldar- 
riaga 2004). Here we show that undetected sources such 
as short-lived quasars or a cluster of faint galaxies, are ex- 
pected to dominate the growth of a galaxy's HII region 
and allow transmission of its Lya photons prior to reion- 
ization. 



2. THE SIZE OF IONIZED REGIONS 

The average mean-free-path for an ionizing photon in a 
fully neutral IGM is much smaller than the size of a cos- 
mological HII region. By balancing the number of ionizing 
photons and neutral atoms we can therefore estimate the 
characteristic physical size of a spherical HII region sur- 
rounding a galaxy with a dark matter halo of mass Afhaio 
(Loeb et al. 2004) 
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where f^, is the star formation efficiency and fcsc is the 
fraction of ionizing photons that escape from the galaxy 
into the surrounding IGM of neutral hydrogen fraction 
xhi. Here we have assumed that the number of ioniz- 
ing photons per baryon incorporated into (Pop-II) stars is 
Y^ - 4000 (Bromm, Kudritzki & Loeb 2001), and that 
Arcion photons are required per baryon for the reioniza- 
tion of the HII region. The value of Aicion depends on the 
number of recombinations and hence on the dumpiness of 
the IGM. The halo mass iVfhaio is a function of its velocity 
dispersion a [sec Eqs. 22-25 in Barkana & Loeb (2001)], 
which we set equal to Vc / \/2 , where Vc is the halo circular 
velocity. 

Similarly, the physical size of the HII region generated 
by an accreting supermassive black-hole of mass Mbh is 
(White et al., 2003) 
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where A^^s is the production rate of ionizing photons dur- 
ing the luminous quasar phase for a black hole of mass 
A^bh = lO^Af©, and /dyn is the fraction of the halo dy- 
namical time during which the quasar shines at its peak 
luminosity. 

There is ample observational data on the number counts 
and clustering of quasars. We therefore compute the size 
of quasar HII regions using a simple, physically-motivated, 
and highly successful model for high redshift quasars that 
has been developed in earlier papers (Wyithe & Loeb 
2003a, 2004). The slow evolution of the quasar clustering 
correlation length with redshift (Groom et al. 2001) implies 
that the black-hole mass scales with halo velocity disper- 
sion independent of redshift. This conclusion is supported 
by direct observations of the Mbh-c relation (Shields et 
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al. 2003). The observed scaling and its redshift depen- 
dence [Mbh oc oc Mi^^f^(l + z)3/2] follow naturally from 
a scenario where black-hole growth is limited by feedback 
over the galaxy dynamical time (Silk & Rees 1998; Wyithe 
& Loeb 2003a). By associating quasar activity with ma- 
jor mergers, it has been shown (Wyithe & Loeb 2003a) 
that when the quasar lifetime equals the galaxy dynamical 
time the observed number counts of quasars arc accurately 
reproduced over a wide range of redshifts, 1.5 ^ z ^ 6.5. 
The size of the quasar generated HIT region (£)q) may 
therefore be evaluated using equation (2) , with a black hole 
mass that is related to the halo velocity dispersion (Wyithe 

& Loeb 2003a) through Mbh = 10^ (c7/54 km s'^f Mq, 
and a lifetime that is governed by a value of /dyn = 0.035 
corresponding to the typical fraction of the virial radius 
occupied by the luminous component of the galaxy. We 
point out that while the association of a total quasar lu- 
minosity with a dark-matter halo must be made via a the- 
oretical model, the model used here correctly predicts the 
number density of high redshift quasars over a wide range 
of luminosities and redshifts. The total number of ionizing 
photons produced by a quasar (with a universal spectrum) 
during its lifetime, reflects a fixed fraction of its black-hole 
rest-mass energy. 

The ratio between the number of ionizing photons pro- 
duced by the central quasar and by stars is 

r =iif " Vrfa-^ (U^y (N*X^ 
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(3xlo4-i) (t^) • 

Since the dependence of Fhaio on halo mass is steeper than 
linear for = const, quasars become increasingly more 
important for higher mass galaxies. Given that both mas- 
sive stars and quasars produce ionizing photons in short- 
lived episodes, the value of Fhaio reflects their net yield 
of ionizing photons. The HII region generated by quasars 
remains as a fossil after the quasar activity ends, since the 
recombination time is longer than the Hubble time at the 
mean IGM density for 2^8. Hence, relic HH regions from 
prior quasar activity must be considered when estimating 
the transmission of Lya photons from high redshift galax- 
ies. In what follows wc show that this remains true even 
when considering the enhanced recombination rate of the 
overdense IGM around massive galaxies. 

3. STROMGREN SPHERES AT HIGH-REDSHIFT 

Prior to reionization, galaxies generate their own 
Stromgren (HII) spheres with characteristic sizes given by 
equations (1) and (2). These spheres expand into a par- 
tially neutral IGM. If the Stromgren sphere is produced 
by a massive galaxy, then the external neutral fraction 
reflects the mean volume filling factor of all small HII re- 
gions produced by its neighboring low-mass galaxies. The 
clustering of galaxies around a bright galaxy could lead to 
overlap of their individual HII regions, even if the rest of 
the IGM is predominantly neutral. As a result the HII re- 
gion surrounding an observed galaxy could be larger than 
expected from equations (1) and (2) due to a cluster of 
faint galaxies that are below the detection threshold of 
the observations. The physical size of the HII region 



surrounding a galaxy of mass Afhaio can be determined self- 
consistently with the mean neutral fraction into which the 
HII region expands, through the relation 

where the sizes of the Stromgren spheres around the cen- 
tral halo of mass Mhaio fmd a galaxy of mass M, respec- 
tively, are Dg^i„ = i^'iLio.q + ^haio,* and Dl^ = ^ + 
D^*- Equation (4) accounts for the higher IGM density 
due to infall around a massive halo through the nonlin- 
ear density contrast AR(i?s, Mhaio) = p/p (where p is the 
mean background density) at a radius R from a halo of 
mass Mhaio (Barkana 2004). The associated density pro- 
file reproduces the appropriate two-point correlation func- 
tion in the linear regime (Barkana 2004) and is evolved 
into the nonlinear regime in spherical shells. The number 
of galaxies surrounding the central halo is found from the 
Press- Schechter (1974) mass function of halos, which in the 
vicinity of a massive halo must be modified relative to the 
background universe. Within the spherical infall model, 
shells with density p at radius R behave as if they arc 
embedded in a closed universe with a matter density pa- 
rameter rim(-R) = p/Pc, where pc = 3{R/R)^/{8ttG) and 
R is the physical radial velocity of the shell (the cosmolog- 
ical constant is negligible at the high redshifts of interest.) 
Thus, dn{R, Mhaio) / dM is the physical number density of 
halos per unit mass computed in a universe with a den- 
sity parameter Q,^(R) (with a normalization and a growth 
factor based on the local shell density). Each HII region 
reaches its final comoving volume when the production 
rate of ionizing photons is balanced by the recombination 
rate interior to that volume. In this phase, the HII volume 
is inversely proportional to the square of the mean density 
times a clumping factor. The clumping is accounted for 
by iVreion, but an additional factor of A^^ needs to be in- 
cluded on the right-hand-side of equation (4) to account for 
recombinations at the mean density within each shell. We 
assume that the most massive galaxy interior to the HII 
region is the one observed. This conservative assumption 
is reflected in the upper limit of the mass integration. The 
value Mmin is the minimum mass (Barkana & Loeb 2001), 
below which gas infall is suppressed from an ionized region 
that has been heated to lO^K. 

Next we specify the free parameters in our model. The 
absorption spectra of the highest redshift quasars indicate 
that the most recent reionization of hydrogen occurred at 
z ~ 6 (Fan ct al. 2001; Wyithe & Loeb 2004). Moreover, 
there is evidence for a temperature rise in the Lya forest 
associated with the reionization of Hell at z ~ 3.5 (The- 
uns et al. 2002). We set the free parameters by requiring 
that hydrogen be rcionizcd by z = 6, and helium be dou- 
bly reionized by z — 3.5. Population-II stars do not have 
spectra that are sufficiently hard to ionize Hell, and so 
in the absence of Population-Ill stars at low redshifts the 
ionization fraction of Hell is determined by quasars. The 
global ionization fractions of intergalactic hydrogen and 
helium are therefore 
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Fig. 1. — The physical size of cosmological HII regions, and the detectabihty of high redshift galaxies. Le/t; The neutral fraction of the 
intergalactic hydrogen vs. redshift [see equation (5)], with the constraints that hydrogen is reionized at 2 = 6 and that helium is doubly- 
reionized at z = 3.5. The gray, left-most curve shows the neutral fraction of Hell. Center. The size of an HII region at (z = 7) due to clustered 
faint sources around a central galaxy with a velocity dispersion a. Also shown (gray curves) are values for the naive estimate of the size of 
the HII region, assuming a fully neutral IGM at mean density [Eq. (1)]. Right: Optical depth for Ly« absorption by the damping wing of 
the distant IGM, at the center of an HII region around a galaxy with a velocity dispersion a. From left to right, the different sets of lines 
correspond to redshifts z = 6.5, 2 = 7 and z = 10, respectively. The different line types mark the four cases A-D (see main text) labeled in 
the middle panel. 
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^ion.HcIII = a^^^-^^M,Ho,q- (5) 

where I?M,Ho,q ~ O.SSZJM.q assuming the median spec- 
tral energy distribution of quasars (Tclfer et al. 2002), 
and that ^ 50% of the photons of sufficient energy to 
double-ionize helium are utilized in the reionization of hy- 
drogen (Wyithe & Loeb 2003b). The ratio N^^5/N,^ion 
may be determined from the requirement that Fion.Hciii 
approach unity at z ~ 3.5. Assuming, for simplicity, that 
-^rcion — coust and that the spectral slope of quasars blue- 
ward of the ionization threshold of hydrogen is specified by 
their median spectral template (Telfer et al. 2002), the ra- 
tio /*/csc/-^roion Can then be determined from the require- 
ment that Fion,Hii approach unity at z ~ 6. Gas clumping 
is expected to increase with cosmic time. By calibrating 
-^rcion bascd On the reionization of helium and hydrogen 
at z ^ 3.5 and 6, we therefore overestimate the signifi- 
cance of recombinations at earlier times and derive con- 
servative results for the minimum transmission of the Lya 
line. We note that even though our model is calibrated 
by requiring that hydrogen reionization be completed at 
z ^ 6 by Pop-II stars and quasars, an early generation of 
massive Pop-Ill stars may be required to produce a first 
partial or full reionization of hydrogen at z 15-20. An 
early reionization would in provide consistency between 
the quasar absorption at z '--^ 6 and the large scale polar- 
ization anisotropics of the microwave background detected 
by the Wilkinson Microwave Anisotropy Probe ( WMAP; 
Wyithe & Loeb 2003b; Cen 2003). 

Reionization histories may now be calculated based on 
equation (5) for different star formation prescriptions, with 
Nj,5 and fi, adjusted to meet the above global require- 
ments. We adopt feedback-regulated prescriptions (Dekel 
& Woo 2002; Wyithe & Loeb 2003a) in which the star for- 
mation efficiency in a galaxy depends on its velocity dis- 
persion Cr, with = for a < (Tmin, f^, = /*,crit (c/cTcrit)^ 
for CTniin < cr < (Tcrit, and = /*,crit for (T > CTcrit • Wc 

consider four cases in this paper: (A) (Jmin = '^a-it — 
10 km s"^ (B) (Tmin = 10 km s-\ a^rit = 120 km s'^ 
(C) cTinin = 0-crit = 50 km s~^; and (D) tTmin = 50 km s~\ 
Ccrit 120 km s~^. The values of (Tniin = 50 km s""'^, 
and (Tniin = 10 km s~^, correspond to the threshold for 
infall of gas from a photo-ionized IGM and the threshold 



for atomic hydrogen cooling, respectively. The value of 
Ccrit = 120 km s~^ is calibrated empirically based on the 
local sample of galaxies (Kauffman et al. 2003; Dekel & 
Woo 2002). We also examine the case of fJcrit = fmin to 
reflect our ignorance about whether a similar feedback op- 
erates at high redshifts, when the dynamical time within 
galaxies becomes comparable to the lifetime of massive 
stars. We have adopted the set of cosmological parame- 
ters determined by WMAP (Spergel et al. 2003), and find 
that a constraint of late reionization at z ~ 6 requires val- 
ues for the parameter combination -/V^./^.^crit/csc/A'rcion of 
11, 2, 22, and 17 for models A-D. 

In the central panel of Figure 1 , we plot the Stromgren 
radius Rs as a function of the velocity dispersion cr of a 
central galaxy at a redshift of z = 7. The four curves cor- 
respond to cases A-D. These calculations of Rs are self- 
consistent with the average neutral fraction of the IGM, 
which is obtained from the reionization histories (Eq. 5), 
and plotted in the left panel of Figure 1. A rapid, stellar 
reionization that completes by z = 6 would imply a neutral 
fraction of xm ~ 0.2-0.5 at z ~ 6.5 — 7, consistent with 
the analysis of the HII regions around the highest redshift 
quasars (Wyithe & Loeb 2004) . At z = 7, the HII regions 
approach a radius of Rs ~ 0.5 — 1 Mpc for a ^ 100 km s~^, 
comparable to the values inferred for recently discovered 
galaxies (Kneib et al. 2004; Rhoads et al. 2004). Galaxies 
with cr ~ 300 km s~^, corresponding to the hosts of the 
luminous z ^ 6 quasars, have i?s ~ 3 — 4Mpc in agree- 
ment with the observed sizes of their Stromgren spheres 
(White et al. 2003). For comparison, we also show by gray 
lines the naive estimate for the size of the HII region using 
only the stars within the central halo, and assuming a fully 
neutral IGM at the mean density (Eq. (1). 

4. SUPPRESSION OF THE LYa EMISSION LINE 

At z ^ 6.5 — 7, the ionized regions around galaxies with 
cr 100 km s~^ are sufficiently large to allow transmis- 
sion of their hya emission line despite the red damping 
wing of hya absorption from the surrounding IGM. It is 
therefore not surprising that these galaxies are observed, 
even if they reside in an IGM that is significantly neutral 
globally. The optical depth Tdamp to Lya absorption at 
the galaxy redshift (Miralda-Escude 1998) is plotted as a 
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Fig. 2. — Comparisons between the contributions to the HII region of stars relative to quasars, and of the central galaxy relative to the 
surrounding cluster of galaxies. Left: The ratio between the number of ionizing photons contributed to the HII region by quasars and by 
stars as a function of the velocity dispersion of the host galaxy [Eq. (7)]. Right: The fraction of the HII region volume contributed by the 
central galaxy as a function of the velocity dispersion of the host galaxy. In each panel the different line types correspond to cases A-D as 
in Figure 1 for galaxies at z = 7. 



function of a in the right hand panel of Figure 1 for galax- 
ies at z ~ 6.5, z ~ 7 and z ~ 10. Here the size of the HII 
region is based on equation (4), while the neutral fraction 
outside the HII region was computed from equation (5). 
At 2; ~ 10, a galaxy must possess a }^ 200 — 300 km s^^ in 
order to have Tdamp < 1- However at z = 6.5 — 7, Tdamp ^ 1 
for galaxies with cr 70—100 km s^^. Hence, observations 
of galaxies at z ~ 6.5-7 (Kodaira et al., 2003; Rhoads et 
al. 2004; Kurk et al. 2004; Stern et al. 2004) are fully 
consistent with reionization completing at z ~ 6. The re- 
quirement that HII and Helll be reionized at z ~ 6 and 
z ~ 3.5 respectively, make this conclusion robust against 
variations in /*,crit or 

5. CONTRIBUTIONS TO THE IONIZING FLUX 

We now explicitly examine the individual contributions 
from quasars and from the clustering of faint sources to 
the HII region around massive galaxies. First, we find the 
ratio between the quasar and stellar fluxes within the HII 
region, Ftot- The flux ratio Fhaio from the central halo is 
obtained from equation (3). The mean flux ratio (Fgai) 
from other galaxies within the HII region of radius Rg is 
then 



(rgal) - {Rs - i'halo) 



-1 



X/ ,.R^dR T'^M^'^^y-^^T, (6) 



and the average flux ratio within the HII region is therefore 

1 tot - . (7) 

The left-hand-panel of Figure 2 shows rtot(o') for cases 
A-D at z = 7. Quasars make a substantial contribution 
to the Stromgren volume around massive galaxies, partic- 
ularly in case-A. However the cluster of smaller galaxies 
surrounding the central halo contributes an overall ioniz- 
ing luminosity that is mostly stellar, resulting in a value 
of Ftot that is smaller than Fhaio- 



The fraction of the HII region excavated by the central 
galaxy is found from the ratio -Dj^aio/-^s- This quantity 
is shown in the right panel of Figure 2 as a function of a 
for galaxies at z = 7. The regions around massive galax- 
ies are predominantly ionized by the surrounding cluster 
of smaller galaxies. The luminous quasars discovered at 
z 6 have, during their active phase, luminosities in ex- 
cess of the time-averaged luminosity of these clusters, and 
are therefore observed as the source of their HII region. 



6. SUMMARY 

The massive galaxies discovered prior to reionization 
should have Stromgren spheres that are much larger than 
naively expected from the UV continuum flux of their 
stars. This increase arises for two reasons: (i) there is an 
additional contribution from earlier quasar activity within 
the central and surrounding galaxies, and (ii) there is a 
biased clustering of faint sources surrounding the central 
galaxy which results in overlap of nearby HII regions at a 
time before overlap is achieved throughout the universe. 
The large extent of the HII regions allow transmission of 
Lya photons from galaxies. We conclude that the dis- 
covery of Lya emission from bright galaxies at z > 6.5 
does not indicate a highly ionized IGM. Nevertheless, the 
strength of the Lya damping increases as the galaxy lu- 
minosity decreases or the redshift increases (see the right- 
hand panel of Fig. 1). A neutral IGM would therefore 
suppress the low-luminosity end of the Lya luminosity 
function of galaxies at high redshifts. 
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